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Abstract 

We investigate J/ip photoproduction asymmetries in the framework of the NRQCD 
factorization approach. It is shown that the color octet contribution leads to large uncer- 
tainties in the predicted asymmetries which rules out the possibility to precisely measure 
the gluon polarization in the nucleon through this final state. For small values of the color 
octet parameters being compatible with J/tp photoproduction data it appears possible that 
a measurement of J/tp asymmetries could provide a new test for the NRQCD factorization 
approach, on one hand, or a measurement of the polarized gluon distribution from low 
inelasticity events {z < 0.7), on the other. 
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1 Introduction 

The study of double spin asymmetries in certain final states is one of the most promising ways to 
measure the polarized gluon distribution function in the nucleon. The favorable approach is to 
study the processes which can be calculated in the framework of perturbative QCD, i.e. for which 
the involved production cross section and subprocess asymmetries can be predicted with small 
uncertainties. Recently, the measurement of J/ip photoproduction asymmetries was suggested 
as a tool for direct access to the polarized gluon density |I], the analysis of asymmetries at the 
subprocess level was based on results obtained in the framework of the color singlet model (CSM) 
1^, ^. However, the CSM can not be considered as a sufficiently reliable model to describe heavy 
quarkonium production and decay processes. The measured cross section of prompt J/ip and 
if)' production at large pt at the Tevatron exceeds the CSM predictions by more than one order 
magnitude 0. In photon-proton collisions, even after taking into account next-to-leading order 
corrections , the CSM does not describe satisfactorily different distributions in J / ip production 
0- 

The factorization approach (FA) based on non-relativistic QCD (NRQCD) represents a more 
reliable framework to study heavy quarkonium production and decay processes Q]. In extension 
of the color singlet model (CSM) 0, the NRQCD FA implies that the heavy quarkonium can 
be produced through color octet quark- ant iquark states, as well. The production of a heavy 
quark-antiquark pair can be calculated perturbatively because the hard scale is determined by 
the mass of heavy quark, whereas the evolution of the ((5Q)-pair to the final hadronic state 
is parameterized by nonperturbative matrix elements, {O^). The relative importance of the 
different intermediate states, both color octet and color singlet, is defined by velocity scaling rules 
ip. The color singlet long distance parameters, being related to the quarkonium wave function 
at the origin, can be calculated in the framework of nonrelativistic models, while the color octet 
parameters have to be extracted from experimental data. Once this is done, predictions for 
other processes can be made. Unfortunately, the present-day theoretical uncertainties are still 
too large, preventing a determination of these parameters with enough accuracy to check the 



universality of the NRQCD FA |10, 11, 12, 13 



In the framework of the color octet model the main contribution to J/i}) photoproduction is 
coming from the color octet states ^5'o^'* and ^Pj(8) WM. Throughout this paper we exploit the 



fact that the corresponding long distance matrix elements can not be determined separately from 
the prompt J /ip production data at the Tevatron; it is only possible to extract their combination. 
Spin effects in heavy quarkonium production, such as J /if) [0, |15[ and T polarization [16 



p]7| , pTSl , pj|] in unpolarized experiments and spin asymmetries in hadroproduction [^, ^ can 
provide a possibility to check the NRQCD factorization approach. As ratios of cross sections, 
spin-dependent observables do not strongly depend neither on the mass of heavy quarks, nor on 
the factorization and renormalization scale. These parameters are crucial in the prediction of 
cross sections and they usually give rise to large uncertainties in the extraction of color octet 
parameters . 



In the present paper we consider double spin asymmetries in J/ip photoproduction in the 
framework of the NRQCD FA. We investigate the uncertainties depending on the long distance 
color octet parameters and the possibility to use J /ip production asymmetries for a determina- 
tion of the polarized gluon density. We will argue that a measurement of the polarized gluon 
density via double spin asymmetries in charmonium photoproduction can potentially provide 
an additional test of the NRQCD FA, or alternatively, the polarized gluon distribution can be 
measured once the FA is confirmed. 
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In the next section we consider J/ip photoproduction and electroproduction asymmetries at 
the subprocess level The uncertainties in the color octet long distance parameters are discussed 
in section 3. The double spin asymmetries for various sets of long distance matrix elements and 
different parameterizations of polarized parton densities are calculated at y/s = 10 GeV which 
can be considered as representative for both the HERMES and the COMPASS experiment Q. 
We also estimate the expected double spin asymmetries in J/ ip electroproduction for a possible 
polarized electron-proton collider with energy i/i = 25 GeV p^ . 



2 Double Spin Asymmetries for J/ij) Subprocesses 
2.1 J/^Ij Photoproduction 

The double-spin asymmetry An for inclusive J/ip photoproduction is defined as 

j/^ ^ da{^ + t ^ J/^ + X)- da{^ + y ^ J/^ + X) ^ EdAa/d^p 

da{^ +p'+^ J/^ + X) + da{^ + t ^ J/^ + X) Eda/d^p ' ^ ' 

where the arrows over 7 and p denote the helicity projection on the direction of the corresponding 
momenta of the initial particles. 

Throughout this paper we use the standard notation of spectroscopy to define the quantum 
numbers of intermediate quark-antiquark states, (2'S'+i)/^(^^'^)^ Here S denotes the spin of quark- 
antiquark state, J is the total angular momentum, the superscripts (1) and (8) correspond to 
color singlet and color octet states respectively. The long distance parameter ((9fg((^'^+^)Lj)) 

describes the transition from the ^^'^"'"^^Lj '^''-state to a hadron H. 

In lowest order of only color octet states contribute to J/ip production, namely through 
the following subprocesses: 

1+9 ^ (2) 

1 + 9 'P^'l (3) 

These heavy quarkonium states are produced at the kinematical endpoint in 2 — >■ 1 subprocesses, 
i.e. at 2; ~ 1, where z = ■ Pp/P'-fPp. The double-spin asymmetries for these subprocesses 
were calculated in Ref. |21|. Uncontrolled uncertainties at ^ ~ 1 arise from higher order 



terms in the NRQCD velocity expansion; in this particular case the higher order terms are 
enhanced by the kinematical factor 1/(1 — 2;) which makes all predictions unreliable close to the 
kinematical boundary [|12|, In the same kinematics, the J/ip is produced via the diffractive 



process, which is an additional source of uncertainty. Another kinematical characteristic of the 
2^1 subprocesses, Eqs. (2)-(3), is the small transverse momentum of the produced J/V'; the 
transverse momentum of the final hadronic state is determined mainly by internal motion of the 
colliding partons and by the recoil of emitted soft gluons in the hadronization phase. 

All these uncertainties are smaller in J/ip production at larger transverse momenta. The 
leading order contribution in this case comes from 2 — 2 subprocesses which are of the order of 



-f + g^J/ij + g, (4) 

l + 9iq) - (5) 

l + 9iq) ^ 's['^+giq), (6) 

l + 9{q) - 'PP+9iq)- (7) 
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The double spin asymmetries for J / ip production in the framework of CSM, whose corresponding 
subprocess is described by Eq. (4), were studied in Refs. We calculated the cross sections 
for the production of a color octet heavy quark-antiquark pair, Eqs. (5)-(7), for different helicity 
states of initial photon and gluon (quark). The expressions for Act are too long to be presented 
in this papeiQ. Our results after averaging the spins of the initial particles are in agreement with 
those used for the calculation of the J/ ijj photoproduction cross section in the NRQCD FA WW\ 




0.4 0.5 0.6 0.7 0.8 0.9 0.5 1 1.5 2 2.5 3 3.5 4 
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Figure 1: Different contributions to double spin asymmetries in quarkonium photoproduction, a) vs. 
z and b) vs. px, all for the GS LO parameterization (set A). The solid line corresponds to the color 
singlet ^Si state, the dashed line to the color octet ^Sq'^ state, the dash-dotted line to the color octet 
^P?'* states and the dotted line to the color octet '^sf^ state. 



The resulting double spin asymmetries for the photoproduction of various heavy quark- 
antiquark states at y/s = 10 GeV are shown in Fig.'s 1 a and b. They are calculated for 
the GS LO parameterization (set A) of the polarized parton distribution functions (PDF) p4[ . 
For the unpolarized parton distribution functions the GRV LO parameterization evaluated at 

= Pr + was used [^. As can be seen, the asymmetries for the ^Sq^^ (dashed lines) and 
^Pf'^ (dash-dotted lines) octet states are rather different from the CSM prediction (solid line). 
Therefore, the predicted asymmetry for J/ijj production in the NRQCD FA is expected to be 
sensitive to long distance color octet parameters. 



2.2 J/i/j Electroproduction 

We also calculated the double spin asymmetries in heavy quarkonium electroproduction. The 
virtuality of the coHiding photon is an additional large scale in the process. Requiring » 
4m^, some theoretical uncertainties connected with the enhancement of higher order corrections 
in the NRQCD expansion can be avoided. O, EBI. At the same time it is possible to reduce 



the contribution from higher twist effects and from J/ip diffractive production ( for a detailed 
discussion see Ref. [p^] ). 

In the framework of the CSM, J/ ip can be produced in electron-proton collisions only in order 
a'^al subprocesses, e + g ^ e+ (cc) + g. The asymmetries for these subprocesses were calculated 

^ The FORTRAN codes can be provided on request 

^ We are grateful to Michael Kramer for providing us with the analytical results for subprocess cross sections. 
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in Ref. [^. However, the main contribution to the total J/ip electroproduction cross section 
comes from the leading order subprocesses a^a^, similar to (2)-(3), when only color octet states 
can be produced with small transverse momenta in the center-of-mass system of virtual photon 
and proton. In this paper we consider the double spin asymmetries only for these subprocesses. 
The next-to-leading order color octet part, O(a^a^), is smaller than the color singlet model 
prediction and does not change the asymmetry calculated in the framework of the CSM in 
Ref. m significantly. 

In the NRQCD factorization approach in leading order QCD J/il) can be produced through 

The expression for the corresponding cross section was 



two color octet states, ^5*0 and ^Pj. 



obtained in Ref. [p^] : 




yS rriciQ'^ + 4ml) 



3Q' + 28ml {Oi^''i'Po)) 
Q"^ + Ami 



(8) 



32mgg^ {Oj^'^ePo)) 



y 



where G{x) is the gluon distribution in the proton, S is the invariant mass of the electron-proton 
system, and the momentum fraction of the struck parton is given by a; = {Q^ + AmD/yS. 
The corresponding polarized cross section has the form: 



A(T{l+p-^l + J/ij + X) = / -p 



[i-yY 



y 



dy 



yS m.ciQ'^ + Ami) 



(9) 



+ 4^2)2 



mt 



where AG(x) is the polarized gluon distribution in the proton. The double spin asymmetry is 
defined by the ratio All = ^cr{lp)/a{lp). 

In the limit ^ the expressions for the cross sections reduce to the convoluted photopro- 



duction cross sections [|T^, ^ with polarized and unpolarized splitting functions correspondingly: 



lim a(l + p 



lim Aa(l 



p 



J/^ + X) 



a 
2^ 
a 
2^ 



dQ^ 
dQ' 



dy 



dy 



y 

i-{i-yy 



■a{j + p ^ j/ij + x), (10) 

■Aa(7 + P^ J/^ + X). (11) 



3 Matrix Elements and Results 



In J/tp photoproduction the dominant contribution comes from the color singlet ^Si and color 



octet states, ^Sq"^ and '^Pf- The production cross section for the '^Si'' octet state is much smaller 



>(8) 



3 c(8) 



than cross sections for these two octet states |T^. The color singlet matrix element is connected 
to the nonrelativistic wave function at the origin and can be evaluated in the framework of 
potential models. It can also be extracted from the leptonic decay width of the J/ip. The 
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corresponding color octet long distance matrix elements can not be calculated in pQCD and are 
usually extracted from experimental data. For the color octet P-wave states they can be reduced 
to one parameter using the NRQCD spin symmetry relation: 



{O^CPj)) = (2J+l)(Of(^Po)). 



:i2) 



After using these relations there are only two matrix elements left which give the main 



contribution to J/ijj photoproduction, {Og{ So)) and (Og ( Pq))- From the available data of 
prompt J/ip production at the Tevatron it is only possible to extract the combination of these 
matrix elements with rather large uncertainties [pGj |T0| , p!3| . Apparently, at Tevatron energies the 



states ^5*0^'' and ^Pj^^ contribute significantly at small transverse momenta of J/ip, px < 5 GeV. 



>(8) 



Consequently, the fitted value of the combination of matrix elements {Og{^So)) and {Og(^Po)), 



defined as [13 



{Ol'\^Po)) 



(13) 



is very sensitive to effects which change the pt spectrum, such as dependence on factorization 
and/or renormalization scales, choice of parameterization of parton distribution functions etc. 
||10|| . The number r is fitted from experimental data. 

To test the sensitivity of double spin asymmetries on the color octet matrix elements we used 
two sets of parameters extracted from the Tevatron data on prompt J/ if) hadropro duct ion [|13 
as shown in Table I: 





LO 


HO 


scaling 


(ore^i)) 


(7.63 ± 0.54) • 10-1 GeV=^ 


(1.30 ± 0.09) GeV^ 




{oi/^{'s,)) 


(3.94 ± 0.63) • 10-3 GeV^ 


(2.73 ±0.45) ■ 10-=^ GeV^ 






(6.52 ± 0.67) ■ 10-2 GeV^ 


(5.72 ± 1.84) ■ 10-3 GeV^ 




r 


3.47 


3.54 





Table I. Values of the long distance matrix elements from LO and HO-improved fits to the CDF data. 



The LO set of the matrix elements corresponds to the leading-order fit of the prompt J/ip 
hadropro duct ion data at the Tevatron and the obtained value for M//'^ is the largest among 
existing fits in the literature. These values of matrix elements lead to an increase of the photo- 
production cross section as z ^ 1 and dramatically overestimate the HI and ZEUS data on J/ ip 
photoproduction at HERA. 

When fitting of HO set, the higher-order corrections were taken into account only approx- 
imately by considering the initial and final state gluon radiation in the framework of a Monte 
Carlo simulation p7| , |TB|]. Including effectively the primordial transverse momentum of partons 
increases the production cross sections of ^Sj^^ and ^pj^-* color octet states at small values of pt- 
The impact of this effect is a significant decrease of the M^/^ combination by about one order 
of magnitude. It is worth noting that this type of HO corrections changes the photoproduction 



cross sections for the color octet states only by about 20% — 30% at HI and ZEUS energies [13 



The HO parameters, presented in Table I, are the smallest ones from the existing fits and are 
compatible with the HI and ZEUS data on the z-distribution of the J /ip photoproduction cross 
section. We consider this set as an example of small color octet parameters. 

In Fig. 2 a) and b) the expected asymmetries in J /if) production, calculated at ^/s = 10 
GeV for the GS LO (set A) parameterization of polarized PDF's, are shown in dependence on z 
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Figure 2: a) z and b) pt dependence of double spin asymmetries of J jil) photoproduction for LO and 
HO sets of color octet matrix elements. The solid curves correspond to the CSM prediction. The dashed 



curves corresponds to the case when [O'lJ'^ i^Po)) = 0, the dash-dotted ones to {O'^' ^ Q- Sq)) = 0. 



and As can be seen from these fig.'s, the uncertainty in the asymmetry of J/ip production 
is large for the LO set of color octet matrix elements because the values of {O^^^ (^Sq)) and 
{O'l^^ {^Pq)) are not fixed separately but only as the combination (13). For the HO set the 



contribution to the cross section of the and '^Pf^ color octet states is smaller, which 

reduces the deviation from the asymmetry predicted by the CSM especially for lower z-values 
{z < 0.7). To make a choice between the different sets of long distance parameters, more 
accurate measurements at energies lower than that for HI and ZEUS are required in the future. 

The COMPASS experiment at CERN § and later 
possibly a presently discussed polarized electron- 
proton collider (EPIC) at = 25 GeV ^ will 
work at higher luminosities and may hence provide 
such measurements. The gluon polarization itself 
at that time has probably been measured in semi- 
inclusive reactions with less uncertainties, for exam- 
ple in open charm production at COMPASS or in 
dijet and direct photon production at RHIC [^. If 
in this situation, the HO set of parameters should 
fit the data of the lower energy experiments, the 
J/ip production asjTiimetries can be used to test 
the NRQCD FA. Alternatively, based upon the HO 
set of parameters and relying on the NRQCD FA 
approach, the polarized gluon distribution can be 
measured with acceptable theoretical uncertainties, 
if only J/ip events with z < 0.7 will be used. 

Using the HO set of color octet parameters in 
Fig. 3 the double spin asymmetries for J/ip photo- 
production are presented at -y/i = 10 GeV for the 
different sets A, B, C of GS NLO parameterizations for polarized PDF's pi]. We used here 



0.6 






LL 


VS=10GeV 




0.5 






0.4 




GS NLO a 


0.3 






0.2 






0.1 




-.-"^ 


GSNLOb 







GSNLOc ■ 










0.5 1 1.5 


2 2.5 3 3.5 A 



Figure 3: The double spin asymmetries for 
J ftp photoproduction at ^/s = 10 GeV us- 
ing HO long distance parameters versus pT- 
The definition of lines is the same as in the 
previous figure. 



7 



different sets of NLO parameterization to show tlie sensitivity of tlie asymmetries to tlie polar- 
ized gluon distribution function, wfiicli is rattier different for tliese tfiree sets. We note tliat, as 
expected, tlie result for the GS LO set A parameterization is practically the same as for the GS 
NLO set B. We show only the dependence of the asymmetries. As was mentioned above, un- 
controlled corrections are expected at large z from higher order terms in the velocity expansion 
because of ignoring the mass difference between the (cc)-pair and the final state hadron. This 
type of uncertainties can be neglected in the da/dpT distribution because the integration over z 
smears the singular region out 
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|]. As can be seen from Fig. 3, when using the HO set of 
color octet parameters the double spin asymmetries for J/ ip production are rather different even 
for the medium size gluon polarization in the nucleon. The deviation from the GSM prediction 
for double spin asymmetries is not large in this case. 

In Fig. 4 the electroproduction asymmetries are 
shown at ^/s = 25 versus photon virtuality, Q^. As 
can be seen from this figure, the difference between 
asymmetries for the states ^Sq"^ and ^Pf^ is large for 
a large gluon polarization. The 'largest' parameter- 
ization, GS NLO (set A), was used to calculate the 
shown electroproduction asymmetries. One should 
keep in mind that in the leading order electropro- 
duction case large contributions from higher twist 
effects and diffractive J/ip production are expected. 
To suppress these contributions it is necessary to re- 
quire large enough » 4m^ where, on the other 
hand, the cross section falls down rapidly and large 
statistical errors are expected p3[. It is a question of 
attainable statistics at EPIG whether it may be pos- 
sible to achieve reasonable restrictions on the values 
of the color octet parameters from J/ip electropro- 
duction asymmetries. At high energies, a/s > 100 
GeV, which correspond to the polarized HERA op- 
tion the total cross section is large but the asym- 
metries for the commonly used parameterizations GS 
or GRSV 1^ |3^ are expected to be very small, less 
than 2%. 
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Figure 4: The Double spin asymmetries 
for J/tp electroproduction at ^/s = 25 GeV 
for GS NLO (set a) parameterization. The 
dashed line corresponds to the asymmetry of 

1 (8) 

Sq state and dash-dotted line to the asym- 
metries of ^Pf"^ states. 



4 Conclusions 

In the present paper double spin asymmetries in J/ ip photo- and electroproduction were studied 
in the NRQGD factorization approach. The asymmetries for the color octet states 5*0 and 
^Pj^^ which give large contributions to the J/ip photoproduction cross section, are significantly 
different from the GSM prediction. The values of the corresponding color octet long distance 
matrix elements can presently not be extracted separately from existing data on J / ip hadropro- 
duction. Two sets of long distance parameters were used to show the sensitivity of the expected 
asymmetries to this uncertainty. For the 'large' set of parameters extracted from GDF data |]13[ 
the asymmetries contain very large uncertainties. For the 'small' set of parameters, when higher 
order QGD corrections are partly taken into account, the predicted asymmetries for the two 
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extreme choices of parameters are rather close to the CSM prediction. The 'small' set of long 
distance parameters is preferable to explain HI and ZEUS data on inelastic J /if) photoproduc- 
tion. The planned high luminosity measurements at COMPASS and possibly at EPIC = 25 
GeV) may provide the possibility to fit the color octet parameters at lower energies with 
better accuracy. If the 'small' set of parameters should remain as the prefered one, the J/ip 
production asymmetries can be used to test the NRQCD FA. Here it is anticipated that the 
gluon polarization will have already been measured with small uncertainties through other semi- 
inclusive modes in polarized electron-proton or proton-proton collisions with small uncertainties. 
On the other hand, if the situation with NRQCD FA is clarified from other experiments, such as 
J / ip and bottomonium polarization measurements, the measurement of double spin asymmetries 
can be used to extract the gluon polarization in the nucleon in LO pQCD using events aX z < 0.7. 
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